Ground-state properties of the two-site Hubbard-Holstein model: an exact solution 
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We revisit the two-site Hubbard-Holstein model by using extended phonon coherent states. The 
nontrivial singlet bipolaron is studied exactly in the whole coupling regime. The ground-state 
(GS) energy and the double occupancy probability are calculated. The linear entropy is exploited 
successfully to quantify bipartite entanglement between electrons and their environment phonons, 
displaying a maximum entanglement of the singlet-bipolaron in strong coupling regime. A dramatic 
drop in the crossover regime is observed in the GS fidelity and its susceptibility. The bipolaron 
properties is also characterized classically by correlation functions. It is found that the crossover 
from a two-site to single-site bipolaron is more abrupt and shifts to a larger electron-phonon coupling 
strength as electron-electron Coulomb repulsion increases. 



I. INTRODUCTION 

It has been demonstrated for more than one decade 
that the polaronic effects observed in High- Tc supercon- 
ducting cuprates [T] and the colossal magnorestive nian- 
ganites [3] are relevant to the electron-phonon (e-ph) cou- 
pling in these systems. The famous Holstein molecular 
crystal model [31], where electrons are coupled with local 
phonons, has then revived in recent years. More recently, 
to include the electron-electron (e-e) Coulomb repulsion 
interactions, the Hubbard-Holstein (HH) model i4j] for 
strongly correlated electron systems has made significant 
progress in understanding many-body aspects. As is well 
known, an exact solution to this model in the thermo- 
dynamic limit is impossible, and many approximate ap- 
proaches are then employed. The two-site HH model, 
where two electrons hop between two adjacent lattice 
sites, can be solved analytically It is not only a pro- 
totype of the HH model but also helpful for the better 
understanding to the polaron and bipolaron behavior in 
infinite lattices. 

Recently, an contemporary alternative characteriza- 
tion of the ground state (GS) properties has been focused 
on quantum information tools in terms of quantum en- 
tanglement d 0, S i, [3 and fidelity 0, [ll [H Q [H, 
[TgI . [l7t . These studies will establish somewhat interest- 
ing understanding from the field of quantum information 
theory to condensed matter physics [H, [l^ . Usually, it 
is hard to calculate these quantities due to the lack of 
knowledge on the exact GS wave function. 

The two-site HH model has been previously addressed 
by means of different methods, such as variational 
method [2l|, perturbation theory [l^l, and numer- 
ical diagonalizationjU (^^j. In spite of these efforts 
there remains poor convergence in the intermediate cou- 
pling regime. A reliable treatment in the whole cou- 
pling regime is still needed. Recently, Berciu derived all 
Greens's function analytically for the two-site HH model 
in terms continued fractions [sj. 



In this work, by using extended bosonic coherent states 
[25I . [2^ . [27I , we develop a new exact technique to deal 
with the two-site HH model. The wave function is pro- 
posed explicitly, by which many quantities can be calcu- 
lated directly. The paper is organized as follows. In Sec 
II we introduce the model and describe the approach. In 
Sec HI we calculate the linear entropy , the GS fidelity 
and its susceptibility to study crossover properties from 
quantum information perspective. The static correlation 
function is also evaluated to analyze the GS properties. 
The main conclusions are given in the last section. 



II. MODEL HAMILTONIAN AND EXACT 
SOLUTION 



The Hamiltonian of the two-site HH model takes the 
form 

i,(T (J 

i i.a 

+ 92UJo^ni^^{bi+s + bl^g) +uJo^blbi, (1) 

'i.cr i 

where i{~ 1 or 2) denotes the label of sites, i + 6 = 2 
for i — 1 and vice versa. Cicr{cl^) is the annihilation (cre- 
ation) operator for the electrons and ni,cr{= c\^Cicy) is the 
corresponding number operator at site i with spin cr, be- 
sides rii = rii^ + rii^. e is the unperturbed site potential. 
t is the usual hopping integral. U and V denote the on- 
site and inter-site Coulomb repulsion between electrons 
respectively, gi and 172 denote the on-site and inter-site 
e-ph coupling parameters. bi{bl) is the annihilation (cre- 
ation) operator for phonons corresponding to interatomic 
vibrations at site i. ujq is the phonon frequency and is 
set unit for convenience in the following. 
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Introducing new phonon operators a = {bi + b2)/V^ 
and d = {bi—b2)/V2, the Hamiltonian ([T]) can be written 
into two independent parts {H = Hd + Ha), 



H, = 



d)d + ^ eUi^ - ^ i{c\cr<^2a + cycler) 

U ^ rii^n^i + Vnin2 

i 

9-{ni-n2)id + d^)-n^gl (2) 



and 



(3) 



where a'^ = a'' + ng+, d = a + i^g+, g+ = {gi + (?2)/\/2, 
and g_ — (gi — g2)/V2. Ha describes a shifted oscillator 
and represents lowering of energy, which is a constant 
motion. Hd represents an effective e-ph coupling system 
which phonons are coupled linearly with the electrons. 

In the effective Hamiltonian Hd, among four differ- 
ent bipolaronic states, such as singlet bipolaronic states, 
singlet Anderson bipolaronic states, singlet bipolaronic 
states (antibonding), and triplet states, the singlet bipo- 
laronic state is the most nontrivial to any approaches. 
We will focus on this state in this paper. 

For three singlet normalized electronic states 



and 



C]^|C]^JO)e, C2|C2jO)e and "^(CilCji — L,j^^L,2p 

the singlet bipolaronic state wave function lip) can be 
expressed as 



cl I C: 



^T)|0)e, 



|(pi)cttC^ 



U 



\0)e 



IT ' 



,t J 



2i 



2T^2l 



|0). 



(4) 



where \f2) and [ip^) correspond to phonon states. 

Inserting it into a schrodinger equation for the effective 
Hamiltonian in Eq. three equations are derived by 
comparing the coefficients of cj|cjj0)e, C2|C2jO)e and 
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)|0). 



[A^A + 2e~4{gl+gl) + U]\^i)-V2t\^3)=E\^,) (5) 



[B^B + 2e-4{gl + gl) + U]\ip2)-V2t\^s)=E\^2) (6) 



and \n)B are Fock states of the new bosonic operators 

Ntr- 

Ifl) =^Cn\n)A 

71=0 
Ntr 



Ec„-^(rft+2.g_)V 



\0)ph (8) 



Ntr 



\ip2) 



n B 



Ntr . 

= J2dn^{d^-2g-y 



n2g_d''-2g^_ 



)ph (9) 



As we know that the vacuum state |0)^^^j is just 
a bosonic coherent-state in d{d^) with an eigenvalue 
2(7(— 2f/_)[2^[2^[23|. So this new basis is overcomplete, 
and actually does not involve any truncation in the Fock 
space of d{d'^), which highlights the present approach. It 
is also clear that many-body correlations for bosons are 
essentially included in extended coherent states ([5]) and 
(HI). As usual, the phonon state \lp^) is expanded in a 
complete basis \n), which is the Fock state of 



Ntr 

n=0 



(10) 



Substituting these phonon states ([5]) , Q , and pUj) into 
Eqs. (I5]),( [B]) and (O and left multiplying state 
s(m|, (to I, respectively, the equations can be written as 
iNtr X 'iNtr matrix, and the m'th row is written as 

[m + 2e - 4(g2 + ^f) + u]c^ - V2t^ A{ni\n) fn = Ec^ 

n 

(11) 

[m + 2e-A{gl+gl) + U]dra-V2t^B{m\n)fr, = Ed^ 

n 

(12) 

[to +2e - 2{gi + 52)' + V]f„^ - V2t{J2 Cn{m\n)A 

n 

+ Ed„(TO|n)s) = ^/™ (13) 



[dU+2e-2{gi+g2f + V]\^s)~V2t{\pi) + \^2))--E\^3) 

(7) 

where we have used two displacement transformation 
= d^ + 2g_ and = d^ - 2g_. Note that the hnear 
term for the phonon operator d{d'^) is removed, and two 
new free bosonic field with operator A{A^) and B{B^) 
appear. In the next step, we naturally choose the ba- 
sis in terms of these new operator, instead of c?(rf^), by 
which the phonon states and \ip2) can be expanded 
in complete basis \n)A and |7t.)b respectively, where \n)A 



where ^(TO|n) = {ni\n 
B(m\n) = {■m\n)A = (-1)"!?, 

p-2g^_ ^min[in,n\f Vm!K!(2g- ) "* + 



mn, with Djnn 



and 



^i=0 {m—i)\{n—i)\i\ 

The eigenvalues E and eigenvectors with coefficients 
; d'a and /„ can be exactly solved by diagonalizing the 
above 3Ntr x 3Ntr matrix numerically. 

To obtain the true exact results, in principle, the trun- 
cated number N^j. should be taken to infinity. Fortu- 
nately, finite terms of the singlet state in Eq. (4) are suffi- 
cient to give very accurate results in the whole parameter 
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FIG. 1: (a) GS energy Eo/t of exact solution (solid line) and 
MLFS transformation (dashed line) vs A for U = 0, 3, 6 and 

10. The red hne represents the twice polaron GS energy, FIG. 2: GS energy Eo/t and the probability of the two-site 
which then separates the polaron and bipolaron regime, (b) bipolarons P2 in the inset vs A for V = 0, 2, 4, 6, 8, 10 by choos- 
The probability of the two-site bipolarons P2 for [/ = 0, 3, 6, ing V = 0.2, t = 2.0, e = 0.8. 
and 10. The other parameters are chosen: V = Q, t = 2.0, 
£ = 0. 



range. It should be noted that in the exact diagonal- 
ization in Fock space of original phonon state d 24l|. a 
considerable large phonon number is needed to give rea- 
sonably good results. We believe that we have exactly 
solved this model numerically. 

To show the effectiveness of the present approach, we 
first calculate the GS energy. Fig. 1(a) presents the GS 
energy Eq jt as a function of the on-site effective coupling 
strength A = g\lt by setting ^2 = conveniently. The 
results for the energy by the variational method based on 
the modified Lang-Firsov transformation with a squeez- 
ing phonon state transformations (MLFS) ^] are also 
list. It is observed that the present results are lower than 
the MLFS results [1^ , especially in the intermediate cou- 
pling regime. Comparing with the Fig. 1(a) in Ref. Q, 
we find that the present results for the GS energy are 
consistent with those from the lowest pole of a Green's 
function. 



III. GROUND STATE PROPERTIES 

A. Crossover from two-site to single-site bipolarons 

As shown in Fig. 1(a) that there is two regimes dis- 
tinctly. In the Ref. [5], the energy vs A curves in two 
regimes are fitted by two functions, and the abrupt drop 
of the first-order derivative signals the crossover regime 
from two-site bipolarons to single-site bipolarons. We 
will propose a quantitative criterion. Because the exact 
wave function in the present technique is explicitly given, 
we can calculate the probability of the system that two 
electrons are in two sites by Eq. (4) directly. 
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P2 = Re{ifz\vz) 



(14) 



FIG. 3: GS energy and the former five excited energies Ejt 
(a) and the corresponding probability P2 (b) against A for 
[/ = 6 by choosing V = 0.2, t = 2.0, e = 0.8. 



The probability of the two-site bipolarons P2 is shown 
in Fig. 1(b). At the weak coupling, the two electrons 
prefer to stay in two sites. As the coupling strength in- 
creases, the two electrons tend to occupy the single site. 
Strictly speaking, at weak (strong) coupling, it is not pure 
two-site (single-site) bipolarons, only dominated two-site 
(single-site) bipolarons. The crossover regime is wide 
for weak on-site Coulomb repulsion U . As U increases, 
the crossover shifts to larger A and becomes more sharp. 
However, the jump of the P2 is not observed, indicating 
unlikely a phase transition. 

The inter-site Coulomb repulsion V may favor the for- 
mation of dominated single-site bipoalrons. To show this 
effect, we calculate the GS energy and the probability of 
the two-site bipolarons for several values of V for fixed 
U, which are displayed in Fig. 2. It is clear that the 
crossover shifts to smaller A with increasing V, and be- 
comes more smooth. 

In our approach, we can also get the excited states and 
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corresponding energies. Both the GS and a few low ex- 
cited state bipolaron energies for U = 6 are presented 
in Fig. [3^a). Interesting, we find that the curves for the 
GS and the 1st excited state energies, the 2nd and 3rd 
excited state energies, 4th and 5th excited state energies 
ahuost merge into a single line in the strong coupling 
regime. The corresponding two eigenstates are almost 
degenerated. Actually they are different. The difference 
is too small to be observed. The merged point shift to 
larger A for higher excited states. To show the bipolaron 
behavior in these states, we also calculate the probability 
of the two-site bipolarons P2 as a function of A, which are 
shown in Fig[3fb). The P2 vs A curves also merge into 
a single line exactly at the same point as the energy vs 
A curves. In the GS and 1st excited state, the value of 
P2 decreases monotonously as A. However, in the other 
high excited states, a non-monotonous behavior of the 
P2 is observed. In a fixed value of A, the probability of 
the two-site bipolarons P2 determines the positive repul- 
sion energy for the double occupancy, the energy gained 
for the two electrons in the same site, and the energy 
gained for the hopping of electrons from one site to other 
site. To increase the energy in higher excited states, the 
competition of these energies contributes the complicated 
behavior shown in FiglSjb). 

B. Linear entropy 

To investigate the crossover from two-site to single-site 
bipolarons in quantum information science, we attempt 
to study quantum entanglement between electrons and 
their surrounding phonons by means of the linear en- 
tropy El psi . |2^, which is an alternative measurement 
to indicate the entanglement of a two-site and single-site 
bipolaron. It is defined as 

Ei = l-Trp^, (15) 

where p — Trph{\ip){'4'\) is the reduced density matrix of 
electrons by taking partial trace over the phonon degrees 
of the freedom. 

The normalized GS wave function of the singlet bipo- 
laronic state in Eq.(4) is described as 

|^) = |^i)|l) + |(p2)|2) + 1^3)13), (16) 

where the phonon states satisfy {(pi\ipi) + ((^2|</'2) + 
i^slfs) = 1- And the quantum states |1), |2) and 
1 3) represent the normalized and orthogonalized singlet 
electronic states c||cjj0)e C2|C2jO)e and "^(cl|C2| — 

c^C2|)|0)e respectively. 

The reduced density matrix p can be derived by taking 
partial trace over the phonon degrees of the freedom p = 
Si j=i Therefore the linear entropy can be 

derived simply as 

El = I - [{ifiilifif + {'P2\(P2f + {fsl'Psf 

+ 2i?e((v^i|^2)' + (^l|V'3>' + (^2|^3>')] (17) 
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FIG. 4: Linear entropy Ei vs A for U = 0, 3, 6, 10 by choosing 
V = 0.2, t = 2.0, e = 0.8. The maximum Ei = 0.5 corre- 
sponds to the single-site bipolaron in strong coupling region. 



As plotted in FigHl the linear entropy Ei increases 
smoothly with the e-ph coupling parameters A for dif- 
ferent e-e Coulomb repulsions U. In general, in weak 
e-ph coupling region, two electrons tend to occupy two 
sites, and form so-called two-site bipolaron, resulting a 
low degree of the quantum entanglement in the original 
coupling strength. For a fixed U, as X increases two elec- 
trons become to tightly interact with the same lattice. So 
the entire charge of the electrons and entire deformations 
of the lattices are restricted on one site, leading to the for- 
mation of a single-site bipolaron. However, the two-site 
bipolaron is known as that each electron just interacts 
with phonons of its own site and then the linear entropy is 
expected to be much smaller than that of the singlet-site 
bipolaron. As shown in FigU] the linear entropy reaches 
its maximum Ei = 0.5 in the strong coupling region, 
which displays that the single-site bipolarons are maxi- 
mally entangled. Obviously the crossover point shifts to 
a large value of the coupling strength X as U increases. 
The similar crossover point has been demonstrated by 
the probability of the two-site bipolarons presented in 
Fig. 1(b). Therefore we can say that in the presence of 
quantum correlations the two-site and singlet-site bipo- 
larons are effectively characterized to a high (low) degree 
of bipartite quantum entanglement. 



C. ground state fidelity 

More recently, another temporary effective quantum 
information tools, i.e., the fidelity has been put forward 
to analyze complicated interacting systems from the per- 
spective of the GS wave functions [l3i[I3|! ^^^^ identify- 
ing polaronic crossover behaviors in bosonic system p^. 
A simple expression of the fidelity F{gi, gi + dgi) is given 
just by the modulus of the overlap of two ground states 
corresponding to two different coupling parameters gi 
and gi -|- 5gi , where gi is a tiny perturbation parame- 
ter [1^. It is expected to signal a peak at the critical 
point. The general GS fidelity is written as 
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D. static correlation function 
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FIG. 5: GS fidelity F and its susceptibiUty S versus A for 
U ^ 0,3, 6, 10 by choosing V = 0.2, t = 2.0, e = 0.8. A 
sharper drop is shown at the crossover for larger U — 10(solid 
line). 



F = mgiMigi + Sgi)) 
= \{'Pi{9i)\'Pi{gi + Sgi)) 
+ {f3i9i)\f3{gi + Sgi))\ 



{^2{gi)\(p2{gi 



Sgi)) 
(18) 



where liplgi)) and {ipigi + <^.9i)) Stve two normalized GS 
corresponding to neighboring Hamiltonian parameters. 
In our calculation Sgi = 10"'' is used. While the fidelity 
susceptibility S{gi) is regarded as a more effective tool 
to detect the singularity at the crossover regime, which 
reads 1301 



^(<?i) 



hm [1 

Sgi — >0 



F{gi,gi + 5g,)]/iSgl) 



(19) 



The GS fidelity F and its susceptibility S* as a functions 
of effective coupling strength A (= / fo^' different e- 
e Coulomb repulsion U are displayed in Fig. [5{a) and 
(b). A narrow dramatic drop is observed in the vicinity 
of the bipolaron transition point as a consequence of the 
dramatic change of the structure of the GS wave function, 
implying that the crossover from the two-site to singlet- 
site bipolaron occurs. As shown in Fig. [^a), the values 
of the GS fidelity are approximated to 1. As is known, 
for the quantum phase transition system the GS wave 
functions from different sides of the level-crossing point 
are almost orthogonal [l3|. However, there is absent of 
the level-crossing and then values of peaks of the GS 
fidelity are around to 1 rather than at the crossover 
point. Further evidence for this crossover is given by 
the susceptibility S, which shows a cusp structure in the 
intermediate coupling regime in Fig. [HKb). Obviously, 
the crossover becomes more abrupt and the critical value 
of the e-ph coupling A where the peak appears becomes 
larger as the Coulomb repulsion U increases from to 10. 
So it is illustrated that the transition behavior from two- 
site to singlet-site bipolaron can be effectively detected 
by a singularity of the GS fidelity and its susceptibility. 



Since the bipolaron crossover transition has been in- 
deed known from the quantum information perspective. 
Naturally, we seek classically to discuss this crossover by 
the means of the static on-site and inter-site correlation 
functions (niui) and (niU2), which reveal the spatial ex- 
tent of lattice deformations induced by electrons respec- 
tively. The Ui denotes the lattice deformations on site 
i produced by the electrons and rii is the number oper- 
ator of the electrons. The GS correlation functions are 
written as 



(niui,2) = {ni[±id + d))/2 - 2gni,2]) 



(20) 



The positive (negative) sign and ni (712) are associated 
with {niuij and (riiM2) respectively. By using the ex- 
act GS wave function of the singlet bipolaronic state ob- 
tained above, the correlation functions Eg. ipO)) can be 
expressed as follows 



N 



[niU2 

N 



m=0 

CjyiCm-l + 0.5/„j/„i_i) 

m— 1 



(21) 



N 



[niUi) 



m=0 

N 

CmCm-l + 0.5/„/m-l) 

m— 1 

J2 V^^{c*^C.m+l+0.5fUrn+l)m 



The functions — (niui) and ~{niU2) against the e- 
ph couphng strength A are plotted in Figs. [6] (a) and 
(b). One can observe that —{niui) and —{niU2) increase 
monotonically with A in weak coupling regime. When A 
exceeds a critical value, — (niui) and —{niU2) show dif- 
ferent behaviors. It is clearly shown that two electrons 
are localized on one site in strong coupling region, result- 
ing the singlet-site bipolaron. It also demonstrates that 
the e-e interaction U affects the magnitude of inter-site 
(on-site) deformations (niU2) (— (niui)) and the transi- 
tion point where — (^1^2) reaches maximum shifts to a 
larger Xa.sU increases. 

The nature of the bipolaron crossover can also be inves- 
tigated by the correlation function ((^i— 7i2)(ui— M2)) /gi- 
FigsID^c) shows that the two-site bipolaron regime of the 
coupling strength A is wider for larger U in the weak 
e-ph coupling regime. This crossover behaviors are con- 
sistent with the above analysis on the GS entanglement 
and fidelity. 
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ent states approach. The GS energy is lower than pre- 
vious results and the crossover regime is revealed clearly 
by the electrons occupancy probability. With the ex- 
act GS wave function, we study the crossover properties 
from quantum information perspective based on the lin- 
ear entropy, the GS fidelity and its susceptibility. The 
single-site bipolarons are maximally entangled and the 
bipolaron crossover transition is more abrupt as e-e in- 
teraction U increases. The similar bipolaron crossover 
behaviors are also observed classically in the static cor- 
relation function. The present study may provide some 
insights into the more complicated Hubbard-Holstein sys- 
tems in a infinite chain. 



FIG. 6: Correlation function versus A (a) — (ni?ii), 
(b)-{niU2), (c) — ((m — n2)(iii -U2))/gi for [/ = 0, 3, 6, 10 by 
choosing V = 0.2, t — 2.0, e = 0.8. A more abrupt crossover 
occurs for larger U = 10 (solid line). 



IV. SUMMARY AND DISCUSSION 

In this present work we have solved exactly the two-site 
Hubbard-Holstein model by the extended phonon coher- 
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